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We study a voltage-controlled version of the superconducting flux qubit [Chiorescu et al, Science 
299, 1869 (2003)] and show that full control of qubit rotations on the entire Bloch sphere can be 
achieved. Circuit graph theory is used to study a setup where voltage sources are attached to the two 
superconducting islands formed between the three Josephson junctions in the flux qubit. Applying 
a voltage allows qubit rotations about the y axis, in addition to pure x and z rotations obtained in 
the absence of applied voltages. The orientation and magnitude of the rotation axis on the Bloch 
sphere can be tuned by the gate voltages, the external magnetic flux, and the ratio a between the 
Josephson energies of the junctions via a flux-tunable junction. We compare the single-qubit control 
in the known regime a < 1 with the unexplored range a > 1 and estimate the decoherence due to 
voltage fluctuations. 



I. INTRODUCTION 

Superconducting (SC) circuits can exhibit a great va- 
riety of quantum mechanical phenomena and are studied 
for their potential as devices for quantum information 
processing. Several different circuit implementations of 
a SC quantum bit (qubit) have been investigated both 
theoretically and experimentally [J, Q • 

A prototype of a SC flux qubit, characterized by a 
working regime in which the Josephson energy dominates 
over the charging energy, Ej Ec, has been theoret- 
ically designed and experimentally realized [1, |3, M, S 
0) H showing quantum superposition and coherent 
evolution of two macroscopic states carrying opposite 
persistent currents that represent the qubit states. The 
flux qubit state is related to a magnetic moment, and 
is thus typically controlled via the application of exter- 
nal magnetic fields which create magnetic flux through 
the loop(s) in the circuit. An advantage of flux qubits 
is their relative insensitivity to charge fluctuations that 
can lead to fast decoherence [13, [HI , while magnetic 
fluctuations are typically more benign. 

A second type of SC qubits, the so called charge qubits 
[3, [S Ei [13 1 operates in the limit in which the 
charge energy dominates, Ec ^ Ej, thus being rela- 
tively insensitive to magnetic fluctuations, while having 
a well defined value of the charge on a SC island, in which 
the presence or absence of an extra Cooper pair deter- 
mines the state of the qubit. The intermediate regime in 
which the Josephson and charge energies are compara- 
ble, Ej w Ec, has been investigated and realized in the 
"quantronium" [isj . Another type of qubit is the Joseph- 
son, or phase, qubit, consisting of a single junction [l9|. 

In this paper, we investigate the possibility of enhanc- 
ing the control of a SC flux qubit via the application of 
electrostatic gates i, [H [IJ. We study the flux qubit 
proposed by Orlando et al. While in [3], the effect 
of any applied voltages was kept low in order to avoid 
charge noise, we explore the possibility of making use of 
the off-set gate charge as an additional control variable. 
We define two device parameters. Assuming for simplic- 




FIG. 1: The voltage-biased SC flux qubit (schematic). The 
circuit consists of a SC ring (blue) with three Josephson junc- 
tions Ji, J2, and J3, threaded by an external magnetic flux 
^x- The Josephson energy of the middle junction J3 differs 
from the other two by a factor of a. A voltage bias Vi is ap- 
plied to each of the two islands formed by the three junctions 
via a capacitor d. 



ity two Josephson junctions to have equal Josephson en- 
ergies {Eji = Ej2 = Ej), the first parameter is given 
by the ratio a — Ej^ / Ej between the Josephson energy 
of the third junction and the remaining two junctions. 
The regime of interest here is 0.5 < a < 1.5 although in 
principle larger values are possible. The second param- 
eter is the ratio between the Josephson energy and the 
charging energy, Ej / Ec which for flux qubits is typically 
about 10 or larger. We analyze the role of these parame- 
ters in detail and, in addition to the well-studied regime 
a < 1, also explore the opposite regime a > 1. Particu- 
lar effort is spent looking for a single-qubit Hamiltonian 
in which an effective pseudo-magnetic field couples to all 
three components of the pseudo-spin represented by the 
circuit. A charge qubit in which a Uy term in the single- 
qubit Hamiltonian has been proposed in (2^ . The possi- 
bility of changing the relative phase of the qubit states, 
together with the capability to flip them, allows full con- 
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trol over the qubit. Full control on the Bloch sphere is 
thought to be very useful in the field of adiabatic quan- 
tum computation [13, [13, 111. 

Circuit theory provides us with a systematic and uni- 
versal method for analyzing any electrical circuit that 
can be represented by lumped elements [H, [S^, HI, HI] ■ 
Through the language of a graph theoretic formalism, 
Kirchhoff's laws and the Hamiltonian of the circuit are 
written in terms of a set of independent canonical coor- 
dinates that can easily be quantized. The formalism of 
[26; 27, 28] is particularly suited for studying circuits con- 
taining superconducting elements, like Josephson junc- 
tions, that are treated as nonlinear inductors. Here, we 
make use of the extended circuit theory that accounts for 
charging effects and can be applied both for charge and 
flux qubits [27I ]. 

Our main result is the identification of the parameter 
range for a and Ej/Ec in the voltage-controlled flux 
qubit in which the single qubit Hamiltonian acquires a ay 
term in addition to the Ux and terms, thus allowing full 
control of the qubit rotations on the Bloch sphere. In this 
regime, we compute the dependence of the single-qubit 
Hamiltonian on the applied voltages Vi and V2. For the 
quantitative analysis of the qubit dynamics we calculate 
the tunneling amplitudes appearing in the Hamiltonian 
as functions of the device parameters. 

The paper is structured as follows. In Section |TT] we 
briefly review circuit theory [H, [23, [H, [2^ and apply it 
to the circuit of Fig.[2]to find its Hamiltonian. Section lllll 
contains the derivation of the effective periodic potential 
in the Born-Oppenheimer approximation. In Section [TVl 
we address the quantum dynamics of the circuit and find 
localized solutions in the periodic potential. In Section Ivl 
we apply Bloch's theory in a tight-binding approximation 
to find general solutions in the presence of a voltage bias. 
Sec. IVII describes the calculation of the tunneling matrix 
elements appearing in the qubit Hamiltonian and their 
dependence on the device parameters a and Ej/Eq- In 
Sec. IVII) we explore the regime {a > 1) and show that 
a full control on the qubit Hamiltonian is feasible. In 
Section IVlIIl we study the decoherence of the qubit due 
to the attached voltage sources. Finally, Sec. IIXI contains 
a summary of our results and conclusions. 



II. THE CIRCUIT 

Here we study a version of the Delft flux qubit [1, 0] 
with an additional voltage control (Fig. [1]). Typically, 
such a qubit circuit also comprises a readout SQUID 
which can be surrounding or attached to the qubit. We 
concentrate on the qubit itself here and do not include 
the SQUID in our analysis because the presence of a 
readout circuit does not alter the analysis and results 
for single-qubit control presented here. A circuit repre- 
sentation of the studied device is shown in Fig. [21 The 
main loop contains three Josephson junctions and the 
loop self-inductance (K), and is threaded by an external 



magnetic flux $2;. The junctions form two SC islands 
to which electrostatic gates with capacitance Ci and C2 
are attached and voltages Vi and V2 are applied. The 
voltage sources represent the new elements in the circuit. 
As long as the junctions are built in such a way that 
the Josephson energy dominates, Ej ^ Ec, the qubit 
is encoded in the orientation of the circulating persistent 
current, as in Rcfs. 3, 7. 

We represent the circuit as the oriented graph Q shown 
in Fig. [2^, consisting of = 8 nodes (black dots) Ui 
[i — 1,...,8) and B — li branches (thin lines) bi 
(i = 1, . . . , 13), in which each branch hi represents one of 
the following lumped circuit elements: a (bare) Joseph- 
son junction J, capacitance C, inductance K, voltage 
source V , and impedance Z . The impedances Zi and Z2 
model the imperfect voltage sources attached from out- 
side to the quantum circuit. Every Josephson junction 
(thick line) consists of 2 branches: a bare Josephson junc- 
tion (J) and the junction capacitance (Cj) as indicated 
in Fig. [2j3. In addition to these two elements, a Joseph- 
son junction can also be combined with a shunt resistance 
[26]. However, these resistances are typically very large 
and can often be neglected; they are not be of interest 
here. The circuit graph Q is divided in two parts. The 
tree is a loop-free subgraph which connects all nodes of 
the circuit and it is represented by solid lines in Fig. [21 
All the branches fi {i = 1, . . . , F) that do not belong to 
the tree are called chords and are represented by dotted 
lines in Fig. [21 In the present case, the number of chords, 
not counting the junction capacitances Cj, is F — 3. 
There can in principle be inductances contained both in 
the tree and in the chords which considerably compli- 
cate the analysis [2^. However, in our case there are no 
inductances in the tree (no L inductances), so that our 
analysis is much simpler than the general one. From now 
on, we make use of the fact that the circuit graph Fig. [21 
has no inductances in its tree. When a chord is added 
to the tree, it gives rise to a unique loop, a fundamental 
loop. In other words, the set of fundamental loops J^i of 
the graph consists of all loops which contain exactly one 
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FIG. 2: a) Circuit of a voltage-biased flux qubit (Fig.[T|. The 
main loop contains three Josephson junctions and a (chord) 
inductance (K). An external magnetic flux threads the SC 
loop. The junctions Ji and J2 are biased by two electrostatic 
gates, representing the main new feature of the circuit. Solid 
lines represent the tree of the circuit graph, while dotted lines 
are the chords, b) Each thick solid line represents a Josephson 
junction shunted by a capacitance Cj. 
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chord fi. The topological information about the graph 
is encoded in the fundamental loop matrix F*-'^-' of the 
circuit {i — 1, . . . , F; j = 1, . . . , B), 



1, if bj G (same direction), 
F\^' = ^ —1, if bj e J^i (opposite direction), (1) 
0, if bj i T,, 



where the direction of the fundamental loop Ti is given 
by the direction of its defining chord ji. The currents I = 
{Ii, . . . , Ib) and the voltages V = {Vi, . . . ,Vb) associated 
with the branches of the graph are divided in into tree 
and chord currents and voltages. 



1= (Itr,l 



ch) 



V= (Vtr,V, 



ch) 



(2) 



With the division into three and chord branches, the fun- 
damental loop matrix assumes the block form 



F(^) = (-F^ 1 1). 



(3) 



We further split up the current and voltage vectors ac- 
cording to the type of branch [27j , 



It,. = (I/,Iy,Iz) 



Ich — (IcjiIcIk) 



Vtr = (Vj,Vy,Vz), Vch = (Vc,,Vc,Vk), 

such that the matrix F acquires the sub-block form. 



1 


Fjc 


FjK 





Fvc 


FvK 





Fzc 


FzK 



(4) 



(5) 



By inspection of Fig. [2l one finds the loop sub-matrices 
of the circuit according to the rule Eq. ([1]), 



' JC 



' vc — 1^ zc 




' JK 



, F 



VK — t> ZK 



(6) 

With Eq. ([3]), Kirchhoff's laws have the compact form 

FIch = -Itr, (7) 

F^Vtr = Vch-*., (8) 

where — {^i, ■ ■ ■ ,^f) is the vector of externally ap- 
plied fiuxes. Only loops with a non-zero inductance are 
susceptible to an external magnetic flux, thus only one 
external flux needs to be considered here, — (0, 0, ^x)- 
The SC phase differences across the junctions ^p — 
{ifi^ ip2, ips) are related to the canonical variables, the 
fluxes through the relation 



(9) 



while the canonically conjugate momenta are the charges 
Q = {Ql^ Q2) on the junction capacitance. Using circuit 




FIG. 3: Plot of the potential U{ip) for tpx = n along the line 
■y^i + <y52 = as a function of — -^(<^i — ^2) for several 
values of a. In the curve for a = 0.5 the two minima are 
degenerate, while for ct > 0.5 they split showing the double 
well. The inset is a density plot of the potential for a = 0.8, 
showing the two minima and the line (/'i + V52 = 0. 



theory [23] and ignoring the dissipative circuit elements 
Zi and Z2 for the moment, we find the following Hamil- 
tonian of the circuit Fig. [21 

Us = i(Q-CyVy)^C-l(Q-CyVy)-fZY(*), 

(10) 

W(*) = -Ejcos27r^ + i*'^Mo* *'^N*x, (H) 

where we have defined cos<^ — (cos (/?!, cos </32, cos (/Sa). 
The Josephson energy matrix is given as 



^ ) =diag(i;j,^j,a^,7j 



(12) 



where $0 = h/2eis the SC quantum of magnetic fiux. We 
assume that the Josephson energies and capacitances of 
the junctions Ji and J2 are equal, Eji = Ej2 = Ej and 
Cji — Cj2 = Cj, and we define the ratio a — Ej^/Ej. 
The capacitance matrices of the circuit are 

Cj = diag(Cj, Cj, Cj3), C = diag(Ci, C2). (13) 

The source voltage vector is defined as Vy = (Vi,V2). 
The derived capacitance matrices C and Cy and the de- 
rived (inverse) inductance matrices Mg and N of Eq. ([1 
are given in the Appendix |^ 



III. BORN-OPPENHEIMER APPROXIMATION 

We consider now the limit in which the chord induc- 
tance K is small compared to the Josephson inductances, 
K <C Lj. By means of the Born-Oppenheimer approxi- 
mation, we derive an effective two-dimensional potential 
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as a function of two "slow" degrees of freedom. Our 
analysis follows closely that of [23|. For K ^ Lj, the 
potential Eq. (fTTjl gives rise to a hard constraint for the 
variables if, in the form of the linear equation 



Mov: + N(^^ 0, 



(14) 



where the external magnetic flux is written as ip^ — 
2tt^x/^o. The general solution of the Eq. (fT4)) . 



<P2 



(15) 



depends on the two variables (pi and only. Thus, in 
the limit of small K, the dynamics is restricted to a plane 
in three-dimensional space. The potential, restricted 
to the plane, is then a function of ipi and Lp2 only 

U{ip) ^ Ej - cos((pi) - cos(v32) - acos((pi - ip2 + Px) 

(16) 

A density plot of U ioi: a — 0.8 as a function of ipi and 
Lp2 is shown in the inset of Fig. [31 The minima of the 
potential are found by solving the equation gradZ^ — 0, 
which yields Q 

simpi = — sin<y92 = — simp*, (17) 

where ip* is the solution of the self-consistent equation 



sinip* — asin(2iy9* + ipx)- 



(18) 



The potential forms two wells whose relative depth is 
determined by the value of the externally applied flux ip^ . 
In order to have a symmetric double well we choose ip^ — 
TT which yields two minima at the points = {ip* , —p*) 
and = (-'/'*)</'*) with ip* = arccos(l/2a) > 0. If 
a > 0.5, then there are two distinct minima. Taking into 
account the periodicity of the potential, a complete set of 
solutions of Eq. is y = ±{(p* , —(p*)'^ +2Tr{n, to), with 
integer n, m. We plot the double well potential between 
the two minima in Fig. [3] for different values of a in the 
symmetric case ip^ = tt. 



IV. QUANTUM DYNAMICS 

In this section, we look for localized solutions of the 
Schrodinger equation — with the Hamiltonian 
of Eq. pil)) . We expand the potential around the two 
minimum configurations, keeping contributions up to the 
second order in tp, and solve the Schrodinger equation in 
these two different points (denoting them L and R for 
left and right). We obtain the quadratic Hamiltonian 



1 r 



Tr -1 



lin; L,R 



(19) 



where the linearized inductance Ljin L r is defined as 



Llin; L,R = Mo 



L/cosv5l,r- 



(20) 



To simplify the kinetic part in Eq. ((T9)) . we perform a 
canonical transformation on the variable # and its con- 
jugate momentum Q [29l |. 



Q = VCQ/y/Z, 



(21) 



where c is an arbitrary unit capacitance (e.g., c = Cj). 
We define the diagonal matrix Jl^ ^ such that it satisfies 



(22) 



where O is an orthogonal matrix that diagonalizes the 
left hand side (Ihs) of Eq. ((22|) . This allows us to further 
simplify the Hamiltonian by making the following canon- 
ical transformation, preserving the Poisson brackets, 

*' = O*, Q' = OQ, (23) 

that leads us to the Hamiltonian, 

HL.R-^(c-^Q'' + *''^f^LR*')- (24) 

In the case of a symmetric potential (when ip^ = tt), the 
matrices Liiii:L,R of the linearized problem are equal. 



Llin;L — Llin;R, and flh — ^Y- 



(25) 



hence we drop the subscript L and R for simplicity. 

We quantize the Hamiltonian by imposing the canoni- 
cal commutation relations. 



Qj] = ihSi^ 



(26) 



where and Qj are the components of the vectors ^ 
and Q respectively. The ground-state wave function is 
the Gaussian, 



^a(<^) 



det 



1/4 



exp 



where a = L, R and 



M = U^] Vco'^noVc. 



(27) 



(28) 



For the wave function overlap integral between the left 
and right state, S = (4'l|^'_r), we find 



S = exp -j -^Aip'^MAifi 



(29) 



where Acp = cp^ — cpi^ = 2 arccos( 1 /2a)(l,— 1) is the dis- 
tance between the right (R) and left (L) potential minima 
(Fig.©. 
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Umax 



Umin 



-2% 27C 

FIG. 4: Density plot of the double well potential U{ipi, ^2) in 
units of Ej for a — 0.8 on a logarithmic scale. The periodicity 
of the potential is evident; each unit cell contains two min- 
ima (black). The primitive vectors of the Bravais lattice are 
denoted ai and a2 while ti and t2 are the tunneling matrix 
elements between the nearest-neighbor minima. 



V. BLOCK THEORY 

Given the periodicity of the problem Eq. pH)) with the 
potential Eq. (flB)) in the Born-Oppenheimer approxima- 
tion, an important question concerns the boundary con- 
ditions of the problem, i.e., the choice of the appropriate 
Hilbert space. The question is whether the domain of cp 
should be the infinite plane or the square T = [—it, tt)^ 
with periodic boundary conditions. This ques tion has 
been discussed extensively in the literature [30l . [Sll . [s^ . 
Since in our case, a shift of ipi or ip2 by 2tt creates a state 
which is physically indistinguishable from the one before 
the shift, we choose the compact domain T and impose 
periodic boundary conditions on the wavefunction. How- 
ever, we are going to extend the domain to the infinite 
domain in order to facilitate the calculation. 



A. The periodic problem 

The approximate solutions constructed in Sec. IIVI are 
a good starting point, but they are insensitive to the 
boundary conditions. However, the boundary conditions 
are essential if finite bias voltages Vy are to be taken into 
account. The problem at hand is defined on the square 
with side 2t: (see inset of Fig. |3]) with periodic boundary 
conditions; i.e., the phases if — {(pi, (^2)^ are in the com- 
pact domain T — [— 7r,7r)^ and the wavefunction at op- 
posite edges needs to be identical, VE'(— tt, 1^2) — ^'(t'', (^2) 
and 5'((/3i,— tt) = ^{(fi,Tr), such that T acquires the 
topology of a torus. If the boundary conditions are ig- 
nored, e.g., in the case where the wavefunction is known 
to be vanishingly small at the boundary, then the bias 



voltages Vy in the Hamiltonian Eq. ((T0|) can be removed 
completely with a gauge transformation and the solutions 
will be independent of Vy . 

We proceed as follows: We first solve the problem 
Eq. (jlOp in the infinite two-dimensional plane and then 
choose those solutions that satisfy the periodic boundary 
conditions and then restrict them to the compact do- 
main T. We choose this approach because the problem 
on the infinite domain is well known: the solutions tpak 
are given by Bloch's theorem for the motion of a particle 
in a crystal and satisfy 

^auiv + 2nm) = e^"™-'' ^,k(v), (30) 

for m = (mi, 7712) with integer mi and m2. The min- 
ima of our potential, Eq. (|16p . define a two-dimensional 
square Bravais lattice with a two-point basis, which looks 
like a sheared hexagonal lattice (although it is a square 
lattice). The lattice and its primitive vectors ai — {2tt, 0) 
and a2 = (0, 2tt) are shown in Fig. [H The lattice basis 
is given by the vectors bL — (0, 0) and bR = 2{ip* ,—ip*). 
Each lattice point can be identified by the Bravais lat- 
tice vector n and the basis index a = L, R. As indi- 
cated above, not all the Bloch functions satisfying the 
Schrodinger equation on the infinite domain have a phys- 
ical meaning, but only those that are also 27r-periodic. In 
the case of zero applied voltage bias, the only value of k 
yielding to a periodic wave function is k = 0. 

B. Tight-binding approximation 

In order to construct approximate Bloch states, we 
first form localized Wannier orbitals (pa by orthonormal- 
izing the locahzed solutions (a = L, R) from Eq. ([77)1 . 
These Wannier orbitals are centered at arbitrary lattice 
points, ipani'^) = — 27711) and satisfy the orthonor- 

mality relations 

(</>an|0/3m) = SapSnrn- (31) 

The Bloch states are then related to the Wannier orbitals 
via a Fourier transform, 

M<P) - E </'-(¥'), (32) 

0an(v) = / dke-2-knV'ak(¥'), (33) 
■/FBZ 

where the integration in Eq. (|33p is over the first Brillouin 
zone (FBZ), i.e., h & [-1/2, 1/2). The label a plays the 
role of the energy band label in Bloch theory. The Bloch 
states ipak form a complete set of orthonormal states in 
k-space, where ki € [—1/2,1/2), 

(V'ak|'0/3q) = Sa(3S{k - q), (34) 
E /"dklV-aOiV-akl =1. (35) 
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For the completeness relation Eq. (|35|) to hold, we must 
sum over all bands a, corresponding to a complete set of 
Wannier functions. Here, in order to describe the low- 
energy physics of the system, we restrict ourselves to the 
two lowest bands a = L, R, related to the left and right 
potential minimum in the unit cell, and neglect higher 
excited states of the double wells. This restriction is jus- 
tified if the energy gap between the lowest two states is 
much smaller than the gap between the two lowest and 
all higher states (see Table We normalize the Bloch 
functions on the unit cell T, 



(36) 



Now we can expand the Hamiltonian in the Bloch func- 
tion basis with Eq. ([55]) . and then apply Eq. ^ 

- ^ J (^kdq|'0ak)('0Qk|H|V'/3q)(V'/3q 



/ dkdq^k^ |V'ak)(V'/3ql, 

a/3 



(37) 



where the approximation in the first line consists in omit- 
ting bands that are energetically higher than a = L, R 
(see above). The matrix elements of the Hamiltonian in 
the Bloch basis are 



(38) 



For fixed k and q, Eq. ([55]) is reduced to a 2 x 2 hermitian 
matrix. The main contributions to Eq. psp stem from 
either tunneling between the two sites in the same unit 
cell (intra-cell) or between site L in one cell and site R 
in an adjacent cell (inter-cell), see Fig. [H For the off- 
diagonal element we can write 



LR 
kq 



E 



— 27r'i(k— q) -n 



bLn\n\cj)nn) 



+ e-2-^«M0L„|H|<^R„-eJ 



\'H\(j)-R.n- 



+62/ 



(39) 



where ei = (1,0) and 62 — (0,1). Due to the lattice 
periodicity, the quantities (see Fig. 2]) 



eo = (0L(fl)n|W|0L(fl)n): 

tl = {4'L{R)n\'H\(t)R(L)n), 

h = (0L(fl)n|W|0fl(L)n-ei> 

= (0L(fl)n|W|0fl(L)n+e2)i 



(40) 
(41) 
(42) 
(43) 



are independent of the lattice site n, and thus from 

Eq. dMl), we find Ti^q - ^(^ - q)Wk'^- We can now 
write the 2x2 Hamiltonian as 



Hk 

A(k) 



If A(k)' 
2 ^(k) 



2 [ti+t2{e 



27Tiki 



(44) 
(45) 



The equality in Eq. ([43]) is due to the invariance of the 
potential under the transformation {ipi, </?2) — *■ —{^2, fi) 
and it is valid also in the (fx case. The eigenvalues 
of the problem are 



^±(k) 



eo±i|A(k)|, 



(46) 



and represent a typical two-band dispersion relation. In 
the case of zero external applied voltage only the k = 
Bloch functions satisfy the correct boundary conditions, 
i.e., are periodic. For k = we recognize the qubit Hamil- 
tonian that, in the symmetric double well case, is given 
by a ax term Q, 



H^eo + {ti+ 2t2) 



C. Effect of a Voltage bias 



(47) 



Now, we study the case with an (nonzero) external bias 
voltage. Given the Bloch function V'ak that satisfies the 
Schrodinger equation for the Hamiltonian Eq. (TTOl) for 
zero applied voltages, = 0, we find for the solution 
wave function for finite voltages Vy 7^ 0, 



— p-iV-Qg/2e 



(48) 



where we have defined the gate charge vector as Qg = 
CyVy. The above statement can be directly verified 
by substituting Uak from Eq. into the Schrodinger 
equation with Eq. (fTO|) while using that ipak solves the 
problem for Vy = 0. The solutions in the presence of an 
applied voltage bias satisfy 

uM<p + 2nn) = e2-"-(k-Q.,/2e) (49) 

For the periodicity of the wave function on the compact 
domain, we have to choose k — Qg/2e. This means that 
Uak is the periodic part of the Bloch function for k ~ 
Qg/2e. By substituting this into Eqs. (gl]) and we 
obtain the qubit Hamiltonian 

H=^ [Re{A)ax + Im(A)a,, + ea,] - • cr, (50) 

where we have also included the effect of a 
(small) bias flux that tilts the double well, 
e ~ 2q!^1 - l/4:a'^Ej{ipx - n), where cr = {<Tx,o-y,az) 
are the Pauli matrices, and 

Re(A) = 2 [ti -I- 2t2COs(7rfc+)cos(7r/c_)] , (51) 
Im(A) = 4t2Cos(7rfc+)sin(7rfc_), (52) 

with k± = (CiVi ± C2V2)/2e. The eigenstates for e = 
are 



|0) = -L(-e-^^|L) + |i?)) 
|1) = -^{e-''\L) + \R)), 



(53) 
(54) 



where tanO — Im(A)/Re(A). In Eq. (|5D|) . we have intro- 
duced the pseudo-field B = (Re(A), Im(A), e). 
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VI. CALCULATION OF h AND fa 

For a quantitative analysis of the singlc-qubit Hamilto- 
nian Eq. (jSOp . we have to calculate the tunneling matrix 
elements ti and t2- In order to do so, we require a set 
of orthonormal Wannier functions on the infinite two- 
dimensional lattice defined by the potential U, Eq. pT|) . 
We start from the non-orthogonal set of Gaussian or- 
bitals l^'an) consisting of the solution Eq. (P7)) . shifted 
by a lattice vector n, 



(55) 



The orthonormalized Wannier functions can be written 
as a linear combination of these Gaussians, 



J2 ^Ml,an|*,d) 



(56) 



To form a complete set of orthonormal functions the fol- 
lowing relation must be satisfied, 



(0a„|0/3m) = (G^Sg) 



an, /3m 



(57) 



where S is the (real and symmetric) overlap matrix, 

/ d<p*a„(<p)^'^m(¥')- 



We solve Eq. ([57]) with 



= g = Vs^. 



(58) 



(59) 



The inverse of S exists due to its positive definiteness. 
The entries of the overlap matrix S are equal to 1 on the 
diagonal, whereas the off-diagonal elements are positive 
and <C 1 because the orbitals 'fan are well localized. We 
define the matrix S''-^-' with all matrix elements ^ 1 via 



5" = l + S-^i 



1 



SlL ShR 



and find, keeping only first order terms in S*'^' 



(60) 



(61) 



Note that 5'll and S'rr have zeros on the diagonal. 

In our tight-binding approximation, we consider 
five unit cells, a center cell with its four near- 
est neighbors, corresponding to the lattice vectors 
{(0,0),(±1,0),(0,±1)}. This means that S and g are 
10 X 10 matrices, which can also be expressed as 2 x 2 
block matrices, each block of dimension 5x5. The 
two largest values are given by si = Rn and S2 = 

n,Rn— ei n,Rn+e2 with the nearest neighbor cell. 

Taking only these two largest overlaps into account, we 
obtain S'll = Sim — and 



S'lr 



/ Si S2 32 \ 

si ' 

si 

S2 si 

V S2 si / 



(62) 



0.75 



0.5 



0.25 



E,/Ec=l 

Ej/Ec=10 

Ej/Ec=15 

Ej/Ec=25 

Ej/Ec=35 



0^ 



0.75 



0.8 



0.85 0.9 

a 



0.95 



FIG. 5: The ratio fa/^i between the tunneling matrix ele- 
ments, plotted as a function of a < 1 for several values of 
Ej/Ec. 



a=l 



0=0.95 




Ej/Ec 



FIG. 6: The ratio t2/ti between the tunneling matrix ele- 
ments, plotted as a function of Ej/Ec for several values of 
Of < 1. 



Having the matrix g and S we can calculate the tun- 
neling matrix 

Taii,f3m = {4>an\'H\4>l3m) = [g^ Tg)a-n,[3rn, (63) 

where the entries of the matrix T are given as 



(64) 



Since both the \'^an) and the \(j)an) states are localized at 
the lattice position n, the matrices T and T both have 
the same non-zero entries as S. The tunneling matrix 
T has the same block form as S with 7ll = ^r = 
eqI and 7lr having the same structure as S'lr with si 
and S2 replaced by ti and t2, given as ti = 7Ln,Rn and 
h = TL„^R„-ei = '?Ln,Rn+e2- Thc Overlaps si and S2, 
together with the transition amplitudes ti and t2, depend 
exponentially on the two parameters a and Ej/Ec- A 
detailed analysis is given below; here, we anticipate the 
approximate relations ti/t2 > 1 if a < 1, ti/t2 < 1 
if a > 1, and ti/t2 ~ 1 if a = 1, and ti/t2 = 1 if 
Ci = C2 = 0. 
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Now, we numerically determine the tunneling matrix 
elements ti and ^2 from Eqs. ([M)) and ([M)) and analyze 
their dependence on the external parameters. This de- 
pendence can then be used to control the qubit Hamil- 
tonian. The external parameters fall into two categories, 
those that can be varied freely, like magnetic fields and 
bias voltages, and the device parameters, that are fixed 
for a specific device. Two main types of device param- 
eters characterize the Hamiltonian: (i) the junction ca- 
pacitance Cj that determines the charging energy Ec — 
e^/2Cj and (ii) the Josephson inductance Lj which de- 
termines the Josephson energy Ej = (<i>o/27r)^/Lj. In 
addition, we have the ratio a = Ej^/Ej. 

The potential U(ip) can be modified in two ways. The 
external magnetic flux = ^^ip^jl-K is responsible for 
the symmetry of the double well within a unit cell and 
can give rise to a CTz term in the single qubit Hamilto- 
nian while a determines the height of the barrier between 
the wells in a cell and between two nearest neighbor unit 
cells. Thus a affects the values of the tunneling ampli- 
tudes between different sites in the lattice. Although a is 
a fixed device parameter for the set-up shown in Fig. [U 
a modified set-up in which the middle junction is made 
flux-tunable has been proposed [l|,[l|; a flux tunable junc- 
tion is achieved by "shunting" the third junction with a 
further junction and using an external magnetic field to 
tune it. 

In the tight-binding picture, the off-diagonal element 
A of the qubit Hamiltonian is a complex quantity that 
depends on the two tunneling amplitudes t\ and whose 
relative strength can be set by a and the ratio Ej/Ec- 
The latter enters as a common factor into the frequencies 
of the Gaussian localized orbitals, determining the size 
of their overlaps and affecting only the energy gap |A|. 
An increase of the value of a implies a decrease of the 
tunneling amplitudes ti and ^2, caused by an increase of 
the height of the barriers. Thus a careful choice of the 
two parameters is crucial in determining the behavior of 
the system. From Eq. ([151) . ih&t if ^2/^1 ^ f 

then A will be (almost) real. In order to obtain a sizable 
imaginary part of A, 12/ ti must be sufficiently large. In 
Fig. O we plot the ratio 12/ ti versus a, for several values 
of the Ej/Ec- Although all the curves approach the 
value ^2/^1 ~ f for a ^ 1, as soon as a < 1, a strong 
variation in 12/ ti is observed for large Ej/Ec- In Fig. [SI 
we plot ^2/^1 versus Ej /Ec for different values of a. For 
a = 1, the curve is almost a constant. In Table [H we 
report a set of quantities calculated by varying both a 
and Ej/Ec, such as to keep the energy gap Aq at zero 
applied voltage of the order of « Q.lEc- 

The parameters of an experimentally realized flux 
qubit (Delft qubit) are a = 0.8 and Ej/Ec = 35 
and are given in the first row of Table [T| In this case, 
the ratio ^2/^1 is very small and the contribution of ^2 is 
negligible. This choice of parameters of the Delft qubit 
therefore does not allow the manifestation of a significant 
CTy term in the single-qubit Hamiltonian, for any value of 
the bias voltage. 




Qi/2e 

FIG. 7: Plot of the real and imaginary part of A as a func- 
tion of Qi/2e = CVi/2e for CV2/2e = 0.5 for a) a = 0.95, 
Ej/Ec = 35; b) a = 0.95, Ej/Ec = 10; and c) a = 1, 
Ej/Ec = 15. 



In Fig. [71 we plot the real and imaginary part of A as a 
function of the applied voltage Vi , expressed in the gate 
charge Qi = CiVi, while keeping the other gate voltage 
fixed such that Q2/2e = C2V2/2e = 0.5. If the real part 
of A can be tuned from a finite value to zero while the 
imaginary part of A remains finite (as in Fig. then 
the pseudo-field B can point along arbitrary angles in the 
equator plane of the Bloch sphere. The magnitude of the 
pseudo-field can be controlled in principle by changing 
a, e.g., with a flux-tunable junction. In Fig. [8l we plot 
the real and imaginary part of A in the case where both 
voltages are varied simultaneously such that Vi = —V2 
as a function of SQ/2e = C{Vi - V2)/2e. In Fig. [i we 
plot the gap |A| as a function of 5Q/2e = C{Vi — V2)/2e 
(solid line) and of {Qi+Q2)/2e = C{Vi+V2)/2e (dashed 
line) for this set of parameters. 



a 


Ej/Ec 


t2/tl 


ti/Ej 
xlO"^ 


t2/Ej 

xlO~^ 


|A|o 
Ej 


|AUi„ 
|A|o 


B12 

|A|o 


0.80 


35 


0.0062 


-2.9 


-1.8 


0.0059 


0.98 


82 


0.85 


30 


0.030 


-1.9 


-5.8 


0.0040 


0.88 


126 


0.90 


25 


0.12 


-1.5 


-18 


0.0037 


0.61 


149 


0.95 


20 


0.39 


-1.5 


-59 


0.0054 


0.12 


116 


1.00 


15 


0.97 


-2.05 


-198 


0.012 





61 


1.05 


10 


1.77 


-4.2 


-740 


0.038 





24 



TABLE L Values of t\, t2, their ratio 12/ ti, the energy gap 
|A|o at zero applied voltage, and the minimum of the gap 
A I mill for a series of values of a and E,j /Ec- In the last col- 
umn we report the ratio of the energy difference E12 between 
the second and first excited state and the qubit gap |A|o. 
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-1 -0.5 0.5 1 



6Q/2e 

FIG. 8: Plot of the real and imaginary part of A as a function 
of 5Q/2e = C{Vi - Vi)l2e for Vi + V2 = choosing a) a = 
0.95, EjIEc = 35, b) Q = 0.95, EjjEc = 10 and c) a = 1, 
EjjEc = 15. 



(Qi+Q2)/2e 

— ' \ 




fll ^ ^ ' ^ ' ^ 

-1 -0.5 0.5 1 

(Qi-Q2)/2e 

FIG. 9: Plot of the gap versus (5Q/2e = 7Cj(Vi - V2)l2e 
(solid line) and 7Cj(Vi + V2)/2e (dashed line), for a. = \ and 
Ej/Ec = 15. In this case both the amplitude of oscillation 
and the cross region of the curves are appreciable. 



finite quantity and can be calculated by minimization of 
equation Eq. (|46p with respect to k. However, for a > 1 
there is a value of the external applied voltage for which 
the gap goes to zero (Fig. [H]). 

We are particularly interested in the regime a > 1. In 
this case ti < t2, i.e., the intra-cell tunneling between two 
minima is inhibited and, with a suitable choice of a, can 
be completely suppressed (Figs. [TUl and [TT|) . In this situ- 
ation, the system can be described by a one-dimensional 
chain in which every even (odd) site is labeled as a "left" 
minimum L while the remaining sites are labeled "right" 
minima R, see Fig. 1121 The tunneling matrix element 
between the sites is ^2 (^i — 0). Note that, due to the 
periodicity of the system, all L (R) sites have to be iden- 
tified with each other, since they describe the same con- 
figuration. 

From Eqs. (jSTjl and ((52|) . we immediately find that, 
for ^1/^2 0, we gain full control of the direction of 
the effective pseudo-field B in the equatorial plane of the 
Bloch sphere, since 

A(fc+,fc_) = 4t2Cos(7rfc+)e"'=-, (65) 

where k± — (CiVi ± C2V2)/2e. The sum and difference 
of the gate charges therefore independently control the 
qubit energy gap and the angle 9 of the pseudo-field, 

|A| = 4|i2C0s(7rfc+)|, e = TTk.. (66) 

VIII. CHARGE DECOHERENCE 

Voltage fiuctuations from imperfect voltage sources 
or other fluctuating charges in the environment lead to 
charge fluctuations on the two islands in the circuit and 
thus to decoherence of the qubit. Moreover, we are con- 
sidering here a situation where the sensitivity to external 
voltages has been deliberately enhanced and therefore it 
can be expected that charge fluctuations cannot be ig- 
nored. An estimate of the decoherence time for the same 



VII. FULL CONTROL FOR a > 1 

The flux qubit realized at Delft |3] operates with a ratio 
a = 0.8 < 1 between the Josephson energies of its junc- 
tions. As shown in Table [H the ratio of tunneling matrix 
elements for this parameter choice is ^2/^1 = 0.0062, thus 
the effect of the applied voltages is negligible. Two other 
regimes for a are interesting, namely a ~ I and a > 1. 

In the former, ti and t2 are approximately equal. In 
this case, can tunnel from a left minimum (L) to a right 
one (R) via both an intra-cell or an inter-cell tunneling 
process with almost equal probability. However, while 
inter-cell tunneling can be controlled via the applied volt- 
ages Vi and V2, allowing superposition with non-zero rel- 
ative phase of the qubit states, the intra-cell transition 
amplitude remains constant, once the parameters a and 
Ej/Ec are fixed, thus leading only to qubit flips. In Ta- 
blelH for each value of a < 1, the minimum of the gap is a 




FIG. 10: The ratio ti/t2 between the tunneling matrix ele- 
ments, plotted as a function of a > 1 for several values of 
Ej/Ec. 
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circuit has been developed in [T^l, where it is found to 
be 0.1s. 

In order to model bias voltage fluctuations, we include 
the two impedances Z\ and (Fig. [5]) in our analysis. 
From circuit theory [l^], we can then obtain a Caldeira- 
Leggett model for the system coupled to its charge envi- 
ronment, 



= 7^5 + Tii 



(67) 



where TYg from Eq. (|T0|) describes the dissipationless el- 
ements of the circuit, and 




2 2 



(68) 



is the Hamiltonian of the degrees of freedom of two in- 
dependent baths of harmonic oscillators that are used to 
model the two impedances, and finally 

UsB^ ^ vcij -Q^CjvXjv, (69) 

j=l,2 V 

describes the system-bath coupling, where mi — 
C-HCi,0)^ and ma = C-i(0,C2)^. The coupling con- 
stants Cjy are related to Zj via the spectral densities 



Jj{uj) = -wReZj(cj) = ^ XI ■ 



-8{lo-uj,,). (70) 



The decoherence rates in the Born-Markov approxima- 
tion are given by [131 



1 

1 



^ J2 |mj- • (0|Q|1)|^ AReZj(A)coth ^ 

i=i,2 



(71) 

|m,-((0|Q|0)-(l|Q|l))|'ReZ,(0)2A:BT. 

(72) 



i=i,2 




Ej/Ec 

FIG. 11: The ratio ti/t2 between the tunneling matrix ele- 
ments, plotted as a function of Ej/Ec for several values of 
a > 1. 



Now we compute the matrix elements of the charge 
operator Q = —2ieV in the |0), |1) basis. Following the 
derivation of the Hamiltonian in Sec. IV B| we start from 

(MQk|Q|"/3k) = -2ek5a0 - 2ie(V'ak|V|V'/3k)- (73) 
The matrix elements of Q between the Bloch states 

(V'aklQIV'/ak) = e'""''"""^ Qom,/3„, (74) 

are given in terms of the matrix elements of V between 
the Wannier functions 



Qan,/3m = -2ei (0q„ | V | (/)/3m) 



/3m ■ 



(75) 



and, in turn, through the ^-matrix, they are expressed 
in terms of the Gaussian states, 

Pan,/3m = (*cn|V|*;3m> ^ A^J^^ ^^S'an,^™ , (76) 



where the matrix M. is defined in Eq. , At^^j^ = 
V/3 ~ ¥'a + 27r(m — n), and the ^-matrix is defined in 

Eq. dsg). 

We only keep the leading matrix elements Si and S2 
in the overlap matrix S when calculating the Q and P 
matrices (see Sec. IVip . Since the largest contributions 
of P are proportional to si and S2, we can use ^ « 1, 
and thus Qan,,3m — Pan, /3m oc So-n./Sm. We cousidcr the 
diagonal term and the off-diagonal term separately and 
obtain. 



(WaklQIWak) 
(wLk|Q|wRk> 



~eiM 



(77) 

siA<^-Hs2(A(^-27rei)e2"*=i 



-|-S2(A</3 -f 27re2) e" 



"27rz/c2 



(78) 




u. 



rmn 



FIG. 12: Density plot of the double well potential U{ipi, 1^2) 
for Q = 1.4, on a logarithmic scale. Two equivalent one- 
dimensional chains with nearest neighbor interaction are high- 
lighted in the figure. 
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where si, S2, Aip = tp^ — ipj^, and the matrix A4 depend 
on a = Ej^/Ej and Ej/Ec- In the qubit basis we find, 

(0|Q|0> - (1|Q|1) = -eAl[sisin(6l)Av. 
+ S2 sin(6' + 2Trfci)(Ay) - 27rei) 
+ S2sin(6l-27rfc2)(A</j + 27re2) , (79) 

(0|Q|1) = ieM\sicos{e)Aip 

+ S2 cos(e' + 27rfci)(Ay> - 27rei) 

+ S2COs(6'-27rfc2)(Ay. + 27re2) , (80) 

where tan^? = ImA/ReA is a function of fci,2 = 
Ci^2Vi,2/2e. Using Eqs. 1^, and HMD we 

can express the decoherence rates in a more explicit way, 



— = 2tt- 



1 

Ti 
1 



2tt 



Ej RcZ 
2kBT ReZ 



^) 52^1(^1,^2), (81) 



T^iVi,V2), (82) 



where S2 , J^i, and are given in the Appendix [Cl J-'i 
and JF^ are periodic functions of the apphed voltages Vi 
and V2 that depend on the parameters a, Ej/Ec, and on 
Si/s2- They can be estimated to be at most of order one, 
depending on the choice of parameters and the applied 
vohages. In Eqs. (HI]) and dHS]) we chose Z w Zi w Z2, 
and Rq = h/e^ is the quantum of resistance. 

In the regime a > 1 we have S2 S> si. For a = 1.4, 
Ej/Ec = 15 and C/Cj 0.02 we find that S2 = 8- lO"*. 
An estimate for T « 100 mK, ReZ « 1 kf] and Ej = 
250 GHz produces decoherence times in the millisecond 
range. 



1 

1 

t1 



Tl{VuV2) 1 

^l,max 6 ms' 
T^{VuV2) 1 



0,inax 



12 ms 



(83) 
(84) 



For some particular values of Vi and V2 the functions Ti 
or vanish, implying that l/Ti ^ or l/T^ — > 0. In 
particular, J^i = for (CiVi, C2V"2)/2e = ±(1/2,0), ± 
(0,1/2), ±(1/4,1/4), ±(1/8,-1/8), ±(3/8, -3/8) in the 
FBZ, and T^^Q for (d^i, C2V"2)/2e ^ (71/2, m/2), ± 
(1/4,-1/4) ± {n,m), with n,m G Z. The two functions 
have a common set of zeros, namely ±(n/2, 0), ±(0, m/2), 
with n,m E Z. In these cases, both 1/Ti, 1/70 0. 

For the regime a < 1 we have that si S> S2 and we 
can neglect terms containing S2- Choosing a = 0.8 and 
Ej/Ec = 35 we find si = 1.3 • 10"^ It follows that the 
decoherence rates are strongly suppressed and an esti- 
mate shows that they are below 1 Hz. This means that 
in this case the main process that causes decoherence is 
not due to the charge degrees of freedom. In fact for the 
Delft qubit 0, that operates in this regime, the dephas- 
ing and the relaxation times caused by other mechanisms 
are much smaller, Ta, = 20 ns and Ti = 900 ns. 



The physical reason for the small decoherence and re- 
laxation rates found here is that, despite the voltage bias, 
we are still dealing with a flux qubit whose states are 
indistinguishable from their charge distribution, as seen 
from Eq. jTZ]). 



IX. RESULTS AND CONCLUSIONS 

By means of circuit theory and a tight-binding approx- 
imation, we have analyzed a voltage-controlled SC flux 
qubit circuit that allows full control of the single-qubit 
Hamiltonian Eq. (j50p . with a^, and terms, in order 
to allow arbitrary single qubit operations. 

One of the main results of this work is the computation 
of the tunneling matrix elements appearing in the single 
qubit Hamiltonian as a function of the device parameters 
a and Ej/Ec- This allowed us to explore new possible 
working regimes of the system, looking for a range of 
parameters for which a full control on qubit rotations is 
feasible. Substantially, the qubit can work in two differ- 
ent regimes, a < 1 and a > 1, showing different features. 
In particular, for a > 1, the pseudo magnetic field B that 
couples to the qubit in the Hamiltonian has a non-zero 
y-component. This allows full control of qubit rotations 
on the Bloch sphere through the applied voltages Vi and 
V2. In fact, in the Hamiltonian, Eq. ([SO)) , the off-diagonal 
term A, given in Eq. (|45| . contains the voltages Vi^2 and 
the sensitivity to Vi_2 is determined by the tunneling pa- 
rameters ti and t2 in Eqs. (gT]), (|i5|) . 

For a < 1, we find ti > <2- The effect of t2, and thus 
of the applied voltages, for the value of parameters of the 
Delft qubit \J], is negligible as shown in Table U but can 
be greatly enhanced for a suitable choice of a and Ej / Ec 
(see Figs. [5]and|ni), thus allowing good control in the real 
and imaginary parts of A, as shown in Eqs. ((5T|) and ((52l) 
and in Figs. [7] and [51 

In the case a > 1, the roles of ti and t2 are inter- 
changed, as shown in Figs. 1101 lll[ and a new regime in 
which a full control of the single-qubit Hamiltonian be- 
comes possible. For a suitable choice of a and Ej/Ec, 
the tunneling parameter ti become vanishingly small, 
giving rise to a simple dependence of A on the voltages, 
as found in Eqs. ([^ and ([S5)) . 

Our analysis is based on the two-level approximation, 
i.e., we assume that we can neglect all high levels be- 
sides the two lowest ones. This approximation is justi- 
fied if the energy gap E12 between the two lowest levels 
and any higher level is sufficiently large, in particular, 
larger than the qubit gap Eqi = |A|. The gap E12 can 
be roughly estimated as the plasma frequency, i.e., the 
smallest of the frequencies of the (anisotropic) harmonic 
oscillator arising from the linearization of the equation 
of motion around the minimum configurations of the po- 
tential. This frequency is given by (also see Appendix [B|) 
WLc = I/VUJL] = y/SEjEc/h. In Table H we report 
the ratio of E12 and the qubit gap |Ao| at zero applied 
voltage. For all parameter values studied, E12 exceeds 
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2| AqI by more than a factor of 20, in many relevant cases 
even by two orders of magnitude, thus justifying the two- 
level approximation. 

Finally, we have studied the decoherence due to charge 
fluctuations of the voltage sources. Our result for the T^^ 
and T^^ rates is given in Eqs. ([8T|) and ([82|) . an estimate 
of which yields a coherence time longer than « 1 ms, 
leading to the conclusion that charge fluctuations are not 
the main source of decoherence, even in the regime in 
which the sensitivity to external voltages is enhanced. 
The coherence of the system is well preserved, since the 
qubit is still essentially a SC flux qubit, i.e., the |0) and 
|1) states have nearly identical charge configurations. 

In conclusion, based on our analysis we find that full 
control of single-qubit operations in a SC flux qubit 
should be feasible, provided that the right choice of the 
device parameters is made. 



via the relation {(pi, ip2, (fis)'^ = P{(pi, (p2)'^ ■ In the case 
of symmetric double well potential, the inductance lin- 
earized matrix LjTjJ.l r is given by 

Lun;L,R = Mo + L7icos<^L,R;i- (B2) 

Because of the symmetry of the potential, we drop the 

subscripts R and L. Applying the matrix V we obtain 
T-l -■pTT-ip 

^lin,P — ' ^lin ' ' 

^rlp-^i r ^""V (B3) 
Lj y — -a a J 
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APPENDIX A: MATRICES C, Cy , Mo, AND N 



In order to simplify the calculation we assume the two 
capacitance Ci and C2 to be equal, Ci = C2 = C and 
define 7 = C/Cj. The projected capacitance matrix 
Cp= V'^CV is then found to be 

y —a 1 + 7 + a J 



The definitions of the derived matrices C, Cy , Mn and 
N that enters the Hamiltonian are given in [26|, [23| for 
the general case. Here we apply the theory and derive the 
matrices for the particular case of the circuit of Fig. [2l 
The derived capacitance matrices are 

C.C,,+ (^°), (Al, 

Cv = (C,0f. (A2) 
The inductance matrices that enter the potential are 

(A3) 



Mo = -FjkF^^, 



N 



K 



(A4) 



and = Mo. For the circuit studied here, we obtain 



Mo = 



1 -1 -1 
-111 
-111 

APPENDIX B: PROJECTED MATRICES 



1 







1-41 











(A5) 



The three-dimensional problem is mapped into a two- 
dimensional one in Sec. IIIII with the matrix 



(Bl) 




In this case, the orthogonal matrices that diagonalize the 
capacitance matrix Cp the linearized inductance matrix 
LiTn.P are identical, Cp = O'^QO and LjT^ip = O'^AO. 
The frequency matrix — dmg{Lj^,uj\\) is given by 



4a2(i+^)2 





40^(1+20+7)^ 



(B5) 



where — 1/LjCj. The matrix M is then diagonal- 
ized by the same orthogonal matrix O and, in the basis 
where it is diagonal, can be written as 



M = 




(B6) 



APPENDIX C: THE FUNCTIONS Ti AND 







We give here an explicit formula for the intra-cell and 
inter-cell overlaps si and S2 as functions of a, Ej/Ec 
and C/Cj, 
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si = exp 



S2 = exp ■ 



E 



J 2 

arccos 



4:V2aEc 



Ej 



^(4a2-l){l + 2a + C/Cj)|, 



16Ec 



1 + C/Cj 



2a 



vr — 2 arccos 



1 

2^ 



(4a2-l)(l + 2a + C/Cj) 
2a 



Through these quantities we can express J^i and J-^/, as functions of fci and fc2, with fc^ = C'iVi/2e, 



Ti{ki,k2) 



|A(fci,fc2)|^^^j^nA(fci,fe) 



2fcBr 



(Cl) 
(C2) 

(C3) 



det^(C) 



^ [7r(Ci,7W22 sin(27rfc2 -9)- C2^Mll sin(27rfci + 9)) 



si 



{C2tMii + CuM22) arccos ( — 1 ( — sin(6l) + sin(27rfci +9)- sin(2fc2 - 9) 



1 2 



(C4) 



where J-f/, is given by J^^, once the sin are replaced by cos. Cij and Mij are the entries of the matrices C and M 
defined in Appendix [b1 The gap |A| and the relative phase between the states |0) and |1) are given by 

|A(fci,/c2)| = 2Y/(ti +2t2C0s[7r(fci - fc2)] cos[7r(fci + /c2)])2 + cos[7r(fci + ^2)] sin[7r(fci - fc2)], (C5) 

2i2COs[7r(A:i +fc2)]sin[7r(fci -/C2)] 

tanty = ;— — 1—7- — -T. (Co) 

<i + 2^2 cos[7r(/ei — k2)\ cos[7r(/ci + /(;2jJ 
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